Three-body non- additive forces between spin-polarized alkali atoms 
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Three-body non-additive forces in systems of three spin-polarized alkali atoms (Li, Na, K, Rb and 
Cs) are investigated using high-level ab initio calculations. The non-additive forces are found to be 
large, especially near the equilateral equilibrium geometries. For Li, they increase the three-atom 
potential well depth by a factor of 4 and reduce the equilibrium interatomic distance by 0.9 A. The 
non-additive forces originate principally from chemical bonding arising from sp mixing effects. 

PACS numbers: 34.20.Mq,31.50.Bc 



There is at present great interest in the properties of 
Bose-Einstein condensates (BECs) formed in dilute gases 
of alkali atoms, and in particular in molecule forma- 
tion in condensates by processes such as photoassocia- 
tion [El @, p|, magnetic tuning of Feshbach resonances 
H and three-body recombination H Once molecules 
have been formed, their fate depends largely on collisional 
processes such as inelastic and reactive scattering, which 
can release kinetic energy and result in the molecules be- 
ing ejected from the trap. Calculations on such processes 
require potential energy surfaces for the three-atom sys- 
tem. In the absence of better information, calculations 
on three-body recombination have mostly used pairwise- 
additive three-atom potentials based on atom-atom pair 
potentials M . The rationale for this is that spin-polarized 
alkali atoms are "honorary rare gas atoms" , and that the 
binding between them is dominated by dispersion forces, 
which are nearly pair wise additive. However, it is known 
that non-additive forces are significant in spin-polarized 
Na3 S, and we have recently shown || that such terms 
can affect ultra-low-energy cross sections for the process 
Na + N& 2 (v = 1) — > Na + Na 2 (u = 0) by at least a factor 
of 10. 

The purpose of the present paper is to investigate the 
magnitude of the nonadditive terms for the complete se- 
ries of homonuclear alkali trimers. We focus on poten- 
tial energy surfaces for quartet states of the three-atom 
systems, corresponding to interaction of spin-polarized 
atoms. These are the surfaces that are most important 
in condensates, and are also the ones for which pairwise 
additivity appears to be a sensible first approximation. 
The corresponding doublet surfaces involve strong chem- 
ical bonding, and are complicated by the presence of con- 
ical intersections at equilateral geometries juj [fl| . 

The approach we have taken is to perform ab initio 
calculations using a single-reference restricted open-shell 
variant [|l2| of the coupled cluster method (ll| with single, 
double and noniterative triple excitations [RCCSD(T)]. 
All the calculations were performed using the MOLPRO 
package [O]. The three-atom interaction potential can 
be decomposed into a sum of additive and non-additive 



contributions, 



■imcr 



^dimcr {Uj ) + V3 (ri2 , r 2 3 , 7-13 ) ■ 

(1) 

The full counterpoise correction of Boys and Bernardi 
|]l5f was employed to compensate for basis set superposi- 
tion error in both dimer and trimer calculations. 

The basis sets used in the ab initio calculations were as 
follows. For Li and Na we used the same basis sets as were 
used by Halls et al. [[l6) and by Gutowski and coworkers 
H [llj respectively. For K, Rb and Cs, we used the small- 
core ECP10MWB, ECP28MWB and ECP46MWB effec- 
tive core potentialsfECPs) of Leininger et al. fig] . These 
quasirelativistic ECPs treat the Is 2 • • -(n-2)s 2 (n- 
2)p 6 (77,-2)d 10 electrons as core and the (n-l)s 2 (n-l)p 6 ns 1 
electrons as valence. To these ECPs, medium-size uncon- 
tracted valence basis sets were added ^2|. The 

resulting atomic electric dipole polarizabilities for Li, K, 
Rb and Cs (165.6, 294.2, 319.34 and 402.38 ag) are in 
excellent agreement with the corresponding experimen- 
tal values H M (164.0 ± 3.4, 292.8 ± 6.1, 319.2 ± 6.1 
and 402.2 ± 8. lag), while that for Na is slightly too high 
(166.3 al compared to 162.7 ± 0.8 ag). 
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FIG. 1: RCCSD(T) interaction energies of spin-polarized al- 
kali dimers. 
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FIG. 2: RCCSD(T) interaction energies of spin-polarized al- 
kali trimers at Dzh geometries (a) total nonadditive poten- 
tials; (b) additive potentials. 



FIG. 3: RCCSD(T) interaction energies of spin-polarized al- 
kali trimers at the Doofc geometries (a) total nonadditive po- 
tentials; (b) additive potentials. 



The RCCSD(T) potential energy curves obtained for 
the a 3 S+ states of the alkali dimers are shown in Fig. 
|. The curve characteristics, V m m = — D e and r c , are 
summarized in Table § The RCCSD(T) results for Li 2 
agree very closely (within 0.1 cm" 1 ) with those obtained 
by Halls et al. using unrestricted quadratic config- 
uration interaction calculations (D c = 334.145 cm" 1 , 
r c = 4.1686 A), which in turn give excellent agreement 
with the RKR curve obtained from optical-optical double 
resonance (OODR) spectra on 7 Li 2 (D e = 333.76 ± 0.02 
cm" 1 , r c = 4.173 A) §| ||. The RCCSD(T) results 
for Na 2 agree very well with those obtained by Gutowski 
using unrestricted coupled-cluster calculations and 
employing the same basis set (D c — 173.926 cm" 1 , 
r c = 5.218 A). Neither our results nor Gutowski's agree 
very well with the RKR curve obtained from OODR spec- 
tra (D e = 175.76±0.35 cm" 1 , r c = 5.108(5) A) @, but 
the accuracy of the RKR curve for Na 2 has been ques- 
tioned Jl|, §§. For K 2 , the RCCSD(T) results are in 
very good agreement with experimental results on 39 K 2 
(D c = 252.74 ±0.12 cm" 1 , r c = 5.7725(20) A) @ 



The RCCSD(T) potentials obtained for the ground 
states of the quartet alkali trimers are shown in Figs. || 
and |^ for equilateral (D^h) and symmetric linear (-Dooh) 
geometries respectively. The full trimer potentials are 
compared with pairwise-additive potentials based on the 
triplet dimer potential in each case. The potential char- 
acteristics, Vaxm = —D c and r c for the global minima 
and Vsu and r sp for the linear saddle points, are listed in 
Table | 

The quartet trimers all have equilibrium interatomic 
distances (at D^h geometries) that are substantially 
shorter than those of the triplet dimers, by an amount 
that decreases steadily down the series from 0.94 A in 
Li3 to 0.59 A in CS3. The three-atom potentials are all 
correspondingly deeper than pairwise sums of dimer po- 
tentials, by a factor that is more than 4 for Li3 but is 1.3 
to 1.5 for the heavier alkalis. The non-additive contri- 
butions V3 to the interaction energies at the equilibrium 
geometries vary from approximately 120% for Li to 50% 
for Cs. These figures are much larger than for systems 
such as the rare gas trimers, where the non-additive con- 
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TABLE I: RCCSD(T) values of r c (A), r sp (A), V mi n = -D e 
(cm -1 ), V sp (cm -1 ), and V3 (cm -1 ) for spin-polarized alkali 
dimers and trimers. 





Dimer 


Trimer D 


ih 


Trimer D 






r c 




r c 


Vini 11 


v 3 


' sp 


V Bp 


v 3 


Li 


4.169 


-334.046 


3.103 


-4022 


-5260 


3.78 


-968 


-354 


Na 


5.214 


-174.025 


4.428 


-837 


-663 


5.10 


-381 


-27 


K 


5.786 


-252.567 


5.084 


-1274 


-831 


5.67 


-569 


-52 


Rb 


6.208 


-221.399 


5.596 


-995 


-513 


6.13 


-483 


-15 


Cs 


6.581 


-246.786 


5.992 


-1139 


-562 


6.52 


-536 


-32 



TABLE II: Natural atomic orbital populations of spin- 
polarized alkali trimers and dimers at the corresponding 
global minima. 



Trimer Dimer 





ns 


np r 


n Pt 


ns 


np z 


Li 


0.743 


0.046 


0.197 


0.992 


0.005 


Na 


0.985 


0.003 


0.009 


0.998 


0.001 


K 


0.949 


0.011 


0.034 


0.995 


0.003 


Rb 


0.975 


0.006 


0.014 


0.996 


0.003 


Cs 


0.947 


0.012 


0.030 


0.995 


0.003 



tributions are closer to 0.5% - 2.5% (31 , 32 and produce 
a weakening rather than a strengthening of the binding. 
However, the figures are quite similar to those for the 
alkaline-earth trimers, where the non-additive contribu- 
tions range from about 100% for Be3 to 60% for Ca 3 Q . 

The RCCSD(T) two-body and three-body interaction 
potentials can be decomposed into self-consistent field 
(SCF) and correlation contributions. The correlation 
contribution dominates at long range, but is overcome 
by the SCF contribution when orbital overlap is signif- 
icant. For triplet alkali dimers, as for rare gas dimers, 
the SCF potentials are repulsive and the main attractive 
forces arise from interatomic correlation (dispersion). 

The qualitative similarity between alkali and rare gas 
atoms does not, however, extend to three-body forces. 
For rare gases, a large part of the three-body energy 
comes from the non-additive dispersion interaction. The 
leading long-range term in this is the Axilrod-Teller- 
Muto (ATM) triple-dipole term ||, ||, which dies off as 
r -3 in each of the interatomic distances; the ATM term 
is repulsive near equilateral configurations but attractive 
near linear configurations. Simulations of rare gas solids 
and liquids using accurate pair potentials and the ATM 
term as the only non-additive contribution have proved 
remarkably accurate [B6j. 

For alkali atoms, by contrast, there is a large attractive 
contribution to the three-body non-additive energy that 
exists even at the SCF level. The SCF values for V3 
at the potential minima are within 20% of the coupled- 
cluster values (~ 110% for Li, - 90% for Na, - 100% 
for K and Rb, and ~ 120% for Cs). This arises because, 
for alkali and alkaline-earth atoms, there are vacant np 
orbitals that lie relatively close to the ns orbitals. The 
np orbitals can form bonding and antibonding molecular 
orbitals (MOs) of the same symmetry (a[ and e') as those 
formed from the ns orbitals. The sets of MOs of the same 
symmetry interact, lowering the energy of the occupied 
MOs and contributing to bonding. In chemical terms, 
this is essentially sp hybridization. 

To verify that this is indeed the mechanism, we have 
carried out natural orbital population analyses |3^ ] of the 
SCF wavefunctions for quartet alkali trimers. The results 



are shown in Table [H]. At geometries, the three 
atoms lie on a circle and the p-orbital populations may 
be separated into radial and tangential parts. For L13, 
the population in radial np r orbitals is 0.046 and that in 
tangential np t orbitals is 0.197. For the other alkalis the 
fractions are much smaller, but still significantly larger 
than the populations of np z orbitals in the corresponding 
triplet dimers. 

We conclude that non-additive terms make substantial 
contributions to the potential energy surfaces for three- 
atom systems involving spin-polarized alkali atoms. For 
Li, the full potential including non-additive terms is a 
factor of 4 deeper than suggested by pairwise additivity. 
For all the alkalis, the non-additive forces are strongest at 
geometries near the equilateral equilibrium configuration. 
A large part of the non-additive forces exists even at the 
SCF level, and arises from sp mixing effects of a type that 
cannot exist in rare gas systems. Non-additive dispersion 
forces are important at long range, but make a relatively 
small contribution around the potential minimum. 

We are currently working on generating complete po- 
tential energy surfaces for the alkali trimers, for use in 
quantum dynamical collision calculations at ultralow en- 
ergies. 
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